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Abstract 

Strongly warped regions, also known as throats, are a common feature of the type IIB 
string theory landscape. If one of the throats is heated during cosmological evolution, the 
energy is subsequently transferred to other throats or to massless fields in the unwarped 
bulk of the Calabi-Yau orientifold. This energy transfer proceeds either by Hawking 
radiation from the black hole horizon in the heated throat or, at later times, by the 
decay of throat-localized Kaluza-Klein states. In both cases, we calculate in a lOd setup 
the energy transfer rate (respectively decay rate) as a function of the AdS scales of the 
throats and of their relative distance. Compared to existing results based on 5d models, 
we find a significant suppression of the energy transfer rates if the size of the embedding 
Calabi-Yau orientifold is much larger than the AdS radii of the throats. This effect can 
be partially compensated by a small distance between the throats. These results are 
relevant, e.g., for the analysis of reheating after brane inflation. Our calculation employs 
the dual gauge theory picture in which each throat is described by a strongly coupled 4d 
gauge theory, the degrees of freedom of which are localized at a certain position in the 
compact space. 



1 Introduction 



Strongly warped regions or throats are a common feature of the landscape of type IIB 
string theory. More specifically, local geometries which are similar to the Klebanov- 
Strassler throat [1] arise naturally in flux compactifications [2] (see also [3]) and the 
distribution of vacua favours geometries with dynamically generated large hierarchies [4]. 
Under certain assumptions, this can even be turned into a prediction for the statistical 
distribution of multi-throat configurations [5]. 

Multi-throat compactifications have been considered earlier [6,7] on the basis of the 
simpler Randall-Sundrum model [8], which realizes the essential features of the Klebanov- 
Strassler throat in a 5d geometry. Cosmological implications of the energy transfer be- 
tween throats have been studied by a number of authors [9-15]. An important motivation 
for the analysis of cosmologies with heated throats comes from the possibility of realizing 
brane inflation in the strongly warped region of the compact manifold [16]. 

In the present paper, we focus on the energy transfer between different throats in 
a given type IIB compactification. If one of the throats is heated during cosmological 
evolution, the energy is subsequently transferred to other throats or to massless fields 
in the unwarped bulk of the Calabi-Yau orientifold. This energy transfer proceeds in 
two ways. If the temperature in a given throat is high enough, it develops a black hole 
horizon [17, 18] and energy is lost by Hawking radiation. When the temperature drops 
below a critical temperature T c , a finite throat undergoes a phase transition during which 
the black hole horizon is replaced by the infrared cutoff region of the throat [18-20]. 
Subsequently, the throat sector contains a non-relativistic gas of Kaluza-Klein (KK) 
modes which decay to other throats in the Calabi-Yau orientifold. 

In both cases, we calculate the energy transfer rate (respectively decay rate) as a 
function of the AdS scales of the throats and of their relative distance. For the decay rate, 
we also demonstrate how to determine its dependence on angular quantum numbers of 
the decaying KK modes. Moreover, we extend the analysis of [6] to a genuine lOd setup 
(for earlier related work see [13,14]). To this end, we consider two AdSsxS 5 throats 
embedded in a 6-dimensional torus. This is a simplified model, but we argue that our 
results remain parametrically correct also for more general geometries. As compared 
to [6], we find a significant suppression of the energy transfer rates if the size of the 
embedding Calabi-Yau orientifold is much larger than the AdS radii of the throats. This 
effect can be partially compensated by a small distance between the throats. These results 
are relevant, e.g., for the analysis of reheating after brane inflation. 

It has been shown in [21,22] that the absorption cross sections for scalars and trans- 
versely polarized gravitons by an AdSsxS 5 throat agree with those by a stack of N 
D3-branes (for appropriate N). The fact that this agreement is exact in spite of the 
use of leading-order perturbation theory in the strongly coupled regime on the gauge 
theory side is explained by a non-renormalization theorem [29]. Motivated by these re- 
sults, we employ the dual gauge theory picture in which each throat is described by 
the world-volume gauge theory on the corresponding stack of D3-branes. The world- 
volume theories on different D3-brane stacks are coupled by the supergravity fields in 
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Figure 1: Two throats with radii R± and R2 separated by a distance A inside a Calabi-Yau 
orientifold of total size L. 



the embedding manifold. The decay and energy transfer rates then follow from the appro- 
priate quantum-field-theory tree-level diagrams. This calculation is considerably simpler 
than the corresponding analysis in the gravity picture, where one has to solve multi- 
dimensional tunneling problems. 

We will use the above equivalence of the gravity and gauge theory picture also for 
non-zero temperature, where the non-renormalization theorem is violated. However, as 
we will show, this only leads to 0(1) uncertainties. The same is true for the generalization 
to the Klebanov-Strassler (approximate AdSsxT 1,1 ) throat. Given that we are anyway 
ignorant about the detailed geometry of the bulk space and of the specific throats which 
may appear in realistic models, we can tolerate this uncertainty. 

We emphasize that, although we refer to throats and the corresponding large- N D- 
brane stacks throughout the text, our results also apply to stacks of fewer branes. This 
may be useful for the analysis of the cosmology of a standard model which resides on 
D-branes in the Calabi-Yau orientifold and heats up the surrounding throats. 

Our paper is organized as follows. In Sect. [2, we derive the energy loss of a heated 
throat to another throat which is separated from the first one by a certain distance A (cf. 
Fig. [1]). This calculation is performed by modelling both throats by stacks of D3-branes 
and replacing the compact space by a torus. It is then straightforward to derive the 
energy transfer rate by summing over the contributions of bulk KK modes coupling to 
both throats. The resulting parametric behaviour ~ I /A 8 of the leading term remains 
valid for more general 6d compact spaces and for more complicated throat geometries. 

Sect. [3] describes an analogous calculation for the decay rate of KK modes localized 
in one throat to fields in a distant throat. In the gauge theory picture, the decaying 
KK modes are represented by glueballs. Thus, we first derive the effective vertex for the 
coupling of these glueballs to bulk fields. After that, the calculation proceeds analogously 
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to that in the previous section. Finally, we compare certain limiting cases of our result 
with calculations in the gravity picture and with formulae from the literature. 

Our conclusions are given in Sect. HJ where we also outline possible applications of 
our results. A relevant integral is evaluated in the Appendix. 



2 Energy transfer between two throats 

Let us consider a compactification manifold containing two throats, one of which is 
heated to a certain temperature T. An interesting quantity for cosmology is the rate of 
energy transfer to the other throat. In the following, we will determine this rate using 
the description of the throats in terms of D-brane stacks. In this picture, a heated throat 
corresponds to a heated world-volume gauge theory. The world-volume theories on the 
two brane stacks are coupled by the supergravity fields in the embedding space. Thus, 
energy transfer between the two throats is, in this picture, due to processes of the type 
shown in Fig. [2J where fields in the thermal plasma on one brane stack scatter into fields 
on the other brane stack. 




Figure 2: Feynman diagram for the scattering of fields on one brane stack into fields on 
another brane stack. 



We will perform the corresponding calculation for a simple example - two semi- 
infinite AdSsxS 5 throats embedded in a 6-dimensional torus of uniform size L. These 
throats are the near-horizon geometries of black 3-branes, which in turn correspond to 
stacks of D3-branes (see e.g. [3,23]). For each throat, the S 5 radius R is related to the 
D-brane number N by 

As it stands, this is not a consistent compactification since negative-charge objects are 
needed to absorb the flux of the branes. However, in the course of our calculation we 
will argue that including these and other objects (e.g. further D-branes) as well as us- 
ing a different embedding manifold and a different throat geometry only leads to 0(1) 
corrections. 

We will restrict our calculation to the mediation by the dilaton, the Ramond-Ramond 
scalar and the graviton polarized parallel to the branes. In the gravity picture these three 
fields satisfy the same wave equation [22]. Correspondingly, in the gauge theory picture 
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their effect in mediating energy transfer is parametrically the same0 Hence, we can 
further restrict our calculation to one of the three fields, which we take to be the dilaton. 
In particular, we will not consider the effect of fermions living in the embedding manifold. 
In fact, in [24] the absorption cross section of dilatinos by 3-branes was calculated and 
found to agree with the result for the dilaton. Therefore, we expect the fermions to give 
parametrically the same contribution as the fields that we consider. 

The (low-energy) world- volume theory on N parallel D3-branes is M = 4 U(iV) super 
Yang-Mills. Its field content is given by the field strength F a p in the adjoint representa- 
tion, six adjoint scalars X % corresponding to the positions of the branes, and fermionic 
superpartners. The coupling between the dilaton and the field strength follows from the 
standard lOd supergravity action with a stack of D3-branes (see e.g. [25]) 
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where, here and below, we work in the lOd Einstein frame. We ignore couplings to 
fermions, since they are proportional to the fermionic equations of motion and thus give 
no contributions to S-matrix elements [22]. Direct couplings between the dilaton and 
the scalars X 1 are absent. Canonically normalizing the dilaton kinetic term and allowing 
for brane fluctuations, we get [21] 
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where (X) is the position of the brane stack. The X % are also defined such that their 
kinetic terms are canonically normalized. As can be seen from Eq. ([3]), couplings involving 

I /2 

the X 1 as well as F a p are suppressed by extra factors of k{ and can therefore be ignored. 



2.1 Energy loss rate to flat lOd space 

Before we proceed, we should check whether a calculation in terms of weakly coupled 
gauge fields is a good approximation in the strongly coupled regime of the gauge theory. 
At zero temperature, this is adequate due to the non-renormalization theorem derived 
in [29]. However, the gauge theory is at finite temperature, which breaks supersymmetry. 
With supersymmetry being broken, the non-renormalization theorem from [29] cannot be 
expected to hold and it is not immediately clear why to trust our calculation. Therefore, 
we analyse a simple example in both the gauge theory and the gravity picture and 
compare the results. Namely, we consider a heated stack of D3-branes in flat lOd space 
which is dual to a non-extremal black 3-brane and calculate the energy loss rate in both 
pictures. 

We model the heated, strongly-coupled gauge theory on the D3-brane stack by a 
thermal plasma of free fields. In principle, one would have to use finite temperature field 

1 This can also be inferred from the relevant part of the DBI action, which couples them to the 
world- volume theories on the D3-branes. 
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theory for the calculation of the energy loss rate. However, as we are only interested in 
the correct order of magnitude, we can perform a zero-temperature calculation using a 
thermal particle distribution in the initial state. Following from Eq. ([3]), the cross section 
for scattering of two gauge bosons into one dilaton is 

a ~ 4, s 3 (4) 

up to 0(1) prefactors, where y/s is the energy of the gauge bosons in the center of 
mass frame. From Eq. (jlj), we can calculate the rate of energy loss per world- volume 
of the branes induced by this scattering process. This is done by thermally averaging 
the product of cross section and lost energy, in analogy to the standard calculations of 
reaction rates in a hot plasma [26,27]: 

P=\J d3 h d 3 k 2 /(wi) f{u 2 ) av {u x + u 2 ) . (5) 

Here 

/M = 4vr3(e^ _ i) ( 6 ) 

is the distribution function for the gauge bosons, v is the relative velocity of the colliding 
particles, and T is the temperature of the heated gauge theory. Inserting Eq. (jlj) into 
Eq. (jSJ), we get the energy loss rate due to scattering of one gauge boson species. To get 
the total energy loss rate, we have to sum over all species and polarizations. In a U(iV) 
gauge theory there are N 2 gauge bosons. Thus, there is an extra factor of 2N 2 coming 
from the summation. Using Eq. ([T]) and neglecting prefactors of order one coming from 
the integration in Eq. (jSJ), we get 

p~i? 8 T 13 , (7) 

where R is the AdS scale of the corresponding black 3-brane. 

Energy loss from the non-extremal black 3-brane is due to Hawking radiation emitted 
by its black hole horizon. The corresponding rate per brane world-volume p is given by 
a generalization of the Hawking formula (see e.g. [25]). If we restrict ourselves to the 
dilaton, we get 

d 9 k vua T {uj) 



where v is the velocity of the emitted particles and T is the Hawking temperature of the 
horizon. The absorption cross section 0^(0;) of a dilaton by a non-extremal black 3-brane 
was calculated in [30]. The result is ct(^) = & (u) f(ui/T), where u is the energy of the 
incident dilaton, / is some function of the dimensionless ratio uj/T, and (Jq(uj) ~ u 3 R 8 is 
the absorption cross section by an extremal black 3-brane with AdS scale R which was 
already determined in [21]. Inserting 0^(0;) and performing the integral, we get 

p~i? 8 T 13 . (9) 

Here we have neglected prefactors of order one which come in particular from the inte- 
gration over f(u/T). 
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Both results for the energy loss rate, Eqs. (|7j) and (Q, agree up to 0(1) factors. 
Accordingly, a weak-coupling calculation in the gauge theory picture gives the right order 
of magnitude. The crucial ingredient is the fact that the absorption cross section <jt of 
a dilaton by a non-extremal black 3-brane differs from the zero-temperature absorption 
cross section cr only by a function of A = uj/T. By gauge/gravity duality, this means that 
the gauge boson-dilaton vertex is corrected by a function of A at non-zero temperature|§ 
Accordingly, the cross section for the process in Fig. [2] that we will calculate assuming 
weak coupling and zero temperature has to be corrected by a function of A. However, 
inserting the corrected cross section into Eq. and performing the integral will just 
give a different 0(1) prefactor, which we ignore anyway. 



2.2 Energy transfer rate to a different throat 

Let us now calculate the cross section for the process in Fig. |2j To this end, we need the 
KK expansion of the dilaton in a 6d torus, 

0(x, (X)) = ^ e 2 ™^/ L (10) 

n&fi 

where L is the size of the torus and the expression is already evaluated at the position (X) 
of one brane stack. The mass of the nth KK mode is = 27r\n\/L. Inserting Eq. ffTOl 
into Eq. ([3]) and using « 10 = M{^ 4 , one sees that the vertex for the nth KK mode in Fig. [2] 
is 

M? L* ■ { j 

Here the energy in the center of mass frame of the gauge bosons is denoted by \/s. 
Let (Xi) and (X 2 ) be the positions of the two brane stacks inside the T 6 . If we denote 
the relative distance of the stacks by A = (X 2 ) — (Xi) and introduce the shorthand 
a = 2ttA/L, the matrix element corresponding to the process in Fig. [2] is given by 

„2 pinct 

M V (12) 

We have ignored prefactors of order one. For phenomenological purposes, we can safely 
assume a/s < L^ 1 . Namely, since the energy s/s of the colliding gauge bosons is deter- 
mined by the temperature T of the heated gauge theory, this corresponds to T < L^ 1 . 
If this were not the case, the gauge theory would heat up the compact manifold and the 
geometrical picture would be lost. Following from ^Js < one has s < for n > 
and the contribution of the energy yfs in the propagator can be neglected for all but the 
zero mode. Thus, Eq. (fl2|) simplifies to 

M - V — + (13) 



2 This is also the case if one takes finite-temperature effects properly into account on the gauge 
theory side. 
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where the prime denotes exclusion of n = in the sum. Since the 4d Planck scale is 
determined by Mf ~ Mf L 6 , the last term in Eq. ffl3|) simply reflects the fact that the 
zero mode interacts with gravitational strength. The sum, which would be UV divergent 
in absence of the exponential factor, is dominated by terms with large n. It can therefore 
be approximated by an integral: 

pina i 

d % n— -. (14) 

n 2 a 4 K J 

The r.h. side of Eq. f fT4|) results from the fact that the exponential function oscillates 
quickly for |n| > a -1 (a = |a|), effectively cutting off the integral Jf| More precisely, we 
evaluate a similar but more general integral, which we will need in Sect. 13.21 in the 
Appendix. Equation ( Tl4|) follows from this integral in a particular limit, which is displayed 
in Eq. (j35j. 

Inserting Eq. ffl4|) into Eq. (1131) . we find 

2 

^~]P + ^' (15) 

where A = \A\. For an order-of- magnitude calculation, we can neglect the interference 
term in |.M| 2 . The cross section for the process in Fig. [2] then reads 

3 

ff ~Mf# + M^ fOT ^ S < L ' 1 - ( 16 ) 

Inserting this cross section into Eq. (jSJ), we get the energy loss rate due to scattering of 
one particle species into another particle species. To get the total energy loss rate, we 
have to sum over all initial and final state species and polarizations. Let us denote with 
Ni and the number of colors of the heated gauge theory and the gauge theory that is 
being heated, respectively. The summation then gives extra factors of 2N 2 and 2iVf and 
we get, again neglecting prefactors of order one coming from the integration in Eq. (J5j), 

q fv — I — ?_T 3 H — —T 9 17) 

P M\l A* + M^L^ ' { ' 

Using Eq. ([1]), this can be written in a slightly more compact form. Denoting by R\ 
and i?2 the AdS scales of the corresponding throats, we arrive at the main result of this 
section: 

•^^| T 13 + ^| T 9_ (lg) 

An apparent limitation of our analysis is the assumption of a simple toroidal geometry 
for the embedding space. This assumption was used to determine the spectrum and the 
couplings of higher KK modes (which determine the first term in Eqs. (fl3l) and (fl8|) ). By 
contrast, the coupling of the zero mode (which determines the second term in Eqs. (|T3l 



3 One can see in particular that the sum in Eq. (|13[) is effectively cut off before the geometry of the 
throats becomes relevant, justifying our flat-space approximation. 
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and (j!8p ). depends only on the size of the embedding manifold and not on its geometry. 
To see the relative importance of the terms more clearly, we rewrite Eq. ffl8|) as 




If the throat-to-throat distance is large, A ~ L, the second term dominates (recall that 
LT < 1) and the precise geometry is irrelevant. By contrast, for small throat separation, 
A <C L(LT) 1//2 , the contribution of the KK modes is dominant. In this case, the precise 
geometry of the embedding manifold may in principle be relevant. However, it is then 
natural to assume that the curvature scale in the region between the throats is smaller 
than I /A. Furthermore, as we have already pointed out above, the sum in Eq. (fl~3l is 
dominated by contributions with \n\ ~ L/A, corresponding to masses ma ~ A' 1 . Such 
modes are only sensitive to the geometry at distance scales A in the vicinity of the two 
throats, which we just argued to be approximately flat. Thus, the order of magnitude of 
our result will remain correct in most relevant cases, even if the overall geometry is very 
different from that of a torus. 

In particular, we see that O-planes and further D-brane stacks will not change our 
result as long as they are not too close to the two throats. Moreover, we can apply our 
result to situations with one Klebanov-Strassler throat and one AdSsxS 5 throat or with 
two Klebanov-Strassler throats as long as the curvature scale of the space in between the 
two throats is not much larger than I /A. 

In order for the calculation in terms of gauge fields to be justified, the temperature 
of the heated throat has to be larger than its IR/ confinement scaleQ One can then easily 
see from the gravity picture that the finite length of the Klebanov-Strassler throats will 
not change the result qualitatively. This is obvious for the heated throat since the black 
hole horizon hides the IR region. For the throat to which the energy is transferred, the 
argument is as follows: In the gravity picture, energy transfer is due to Hawking radiation, 
which is emitted by the heated throat and subsequently absorbed by the other throat. 
But only the geometry in the UV region of the throat is important for the absorption 
by (or, equivalently, the tunneling into) that throat. 

3 Decay of KK modes between two throats 

Another interesting quantity for cosmology is the rate with which KK modes localized 
in one throat decay to a different throat. This question has already received significant 
attention in the literature (see [6,9-15]), mainly in the context of reheating after brane- 
antibrane inflation. However, in all cases the calculations were done in the gravity picture, 
whereas we will again (mainly) exploit the gauge theory point of view. This will allow us 
to incorporate easily the dependence on the throat radii and the distance between the 
throats. We compare the results from the literature with ours in Sect. 13.31 

4 Otherwise, the heated throat sector contains a non-relativistic gas of KK modes, whose decay rate 
to the other throat will be determined in Sect. [3] 
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Figure 3: Feynman diagram for the Figure 4: Potential in the effective Schrodinger 
decay of a glueball in one throat equation for the dilaton in a throat, 
into fields in another throat. 

3.1 The glueball decay vertex 

We want to calculate the decay rate of glueballs on one brane stack into two gauge 
fields on another brane stack. As in Sect. El we perform the calculation for two D3-brane 
stacks in a 6-dimensional torus of uniform size L. As before, we can argue that our result 
provides the right order of magnitude also for more general geometries. The Feynman 
diagram for the process is shown in Fig. [3] Due to the non-renormalization theorem 
described in the introduction, we do not have to care whether the decay products will 
arrange into one or more glueballs. The vertex for this part of the diagram is simply 
the one already derived in Eq. However, the other vertex between a dilaton and a 
glueball can not so easily be read off from the Lagrangian. Therefore, we make use of 
the gravity picture to calculate the decay rate in a simpler situation. From this we will 
determine the vertex by demanding that this decay rate agree with the gauge theory 
picture. 

Namely, we consider a dilaton localized in a single AdSs x S 5 throat which is embedded 
into flat lOd space. This is the geometry of an extremal black 3-brane, the metric being 
given by [31] 

ds 2 = fir)- 1 / 2 {-dt 2 + dx\ + dx\ + dx 2 3 ) + f{r) 1 ' 2 (dr 2 + r 2 d^) 

with f(r) = 1 + ^. (20) 

In the near-horizon region, for r <C R, the warp factor reduces to /(r) ~ R 4 /r 4 and 
the geometry is asymptotically AdSsxS 5 . Far from the horizon, for r 3> R, the warp 
factor is f(r) ~ 1 and the geometry is asymptotically flat lOd space. In this situation, 
the AdS/CFT conjecture is based on taking the near- horizon limit r — > and a' — > 0, 
while keeping r/a! fixed. This effectively reduces the geometry to the AdSsxS 5 part. In 
the equivalent description by a stack of D3-branes in flat space, interactions between 
supergravity and the world- volume theory vanish in this limit. This can be seen, e.g., 
from Eq. (|3J), since a' — > implies k 10 ~ g s a' 2 — > at finite g s . One then identifies states 
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in the world- volume gauge theory with eigenmodes of supergravity on AdSsxS 5 . 



In our considerations, however, we want to retain the asymptotically flat part of 
the geometry. What were previously the eigenmodes on AdSsxS 5 will then become part 
of the spectrum of excitations in the full geometry, including excitations outside of the 
throat region. This reflects the fact that, since we do not set a' to zero, the gauge theory 
will interact with the supergravity fields in the embedding space. 

Such nonvanishing interactions lead to the decay of gauge theory states, to which we 
refer as glueballs, into supergravity fields. This glueball decay has a simple counterpart 
on the gravity side. Namely, excitations in the gauge theory correspond to excitations in 
the throat region. The state dual to the glueball will therefore be a wave packet which is 
localized in the throat. Due to the different time evolution of its constituent modes, this 
wave packet will decohere after a certain time (see [15]). Hence, excitations will show up 
in the asymptotically flat region as well, which is the analogue of glueball decay. 

We will now determine the decay rate of a dilaton localized in the throat into flat lOd 
space. To this end, we will assume the throat to be sharply cut off somewhere in the IR. 
Such an AdS 5 x S 5 throat with an IR cutoff might not exist as a solution to supergravity, 
but it can serve as a simple toy model capturing the relevant information. Later on 
we will show how to extend our results to realistic finite throats such as the Klebanov- 
Strassler throat [1]. On the gauge theory side, the cutoff corresponds to a deformation by 
a relevant operator, in which case the gauge theory has a discrete set of glueball states. 

The wave equation for the dilaton is just the Laplace equation in the background 
geometry: 

d M (Vgg MN d N <f>) = o. (21) 

Using Eq. (|20|) . one gets 



dr dr 



+ m z + 



m 2 R 4 1(1 + 4) 



(r) = 0, 



(22) 



where 1(1 + 4) is the eigenvalue of the Laplacian on S 5 and m 2 is the eigenvalue of the 
4d d'Alembertian. We will call m the mass of the excitation. Choosing a new radial 
coordinate z = R 2 /r and introducing a redefined field <j) = z~ 3 ^ 2 (j), we arrive at 



dz 2 



m 



15/4 + /(/ + 4) m 2 R 4 



0. 



(23) 



This has the form of a Schrodinger equation, the potential being given by the term 
in brackets. A schematic plot of this potential is shown in Fig. HI As one can see, a wave 
coming from the near-horizon region (z — * oo corresponding to r — > 0) has to tunnel 
through an effective barrier to reach the asymptotically flat region (z — > corresponding 
to r — ► oo)E| The tunneling probability V has been calculated in [21] (see also [14]) for 

5 As we will see in Sect. 13.31 by using cartesian coordinates for the torus, one again gets a Schrodinger- 
like equation. However, in this case there is no barrier an incoming wave would have to tunnel through. 
Instead, the reflection of a large part of the incoming wave is due to the steepness of the potential well. 
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masses m <C R , 

V ~ (mi?) 8+4 ' , (24) 

where we again neglect prefactors of order one. 

Although this result has been derived for a throat which is infinite in the IR direction, 
we can still use it for a finite throat, as long as the mass m of the wave is not too small. 
For a throat which is cut off in the IR at z = z™, masses are quantized in units of 
m IR = z" 1 such that m n ~ n m m with n an integer The result of Eq. (!24j) can then be 
trusted as long as the wave function is not completely dominated by the unknown IR 
cutoff region. This will be the case if n is sufficiently larger than 1. 

The wave packet describing the glueball can be decomposed into a set of modes 
moving in the IR direction and in the UV direction. If the barrier on the UV side were 
impenetrable, the modes would be reflected entirely on the UV and IR side. However, 
since a small fraction of the incoming flux is able to penetrate the barrier, the wave 
leaks out of the throat. The incoming and outgoing fluxes at the barrier, ji n and j ut? 
determine the tunneling probability V and the decay rate V: 

V = l-j out /j in , T = Jm V. (25) 

Thus, a wave packet localized in the throat will decohere. 

To determine T, we need solutions to Eq. (|2"31) describing waves which are reflected 
back and forth between the UV barrier and the IR end of the throat. From these we can 
calculate the incoming flux j in . We restrict ourselves to the case m <C R~ x . In particular, 
this means that z m ^> R, where z = R corresponds to the beginning of the throat region 
(cf. Eq. (120]) ). For z m _1 ^> R, we can neglect the last two terms in the potential, 
keeping only the constant term m 2 . In this limit, the solution is simply given by plane 
waves: 

(j) ~ A cos mz + B sin mz. (26) 

The approximation is valid for z m > z ^> m" 1 ~ z m /n. If n is not too small, the mode is 
well approximated by a plane wave in a large portion of the throat. Deviations from this 
form for z < z m /n are due to reflection at and tunneling through the effective barrier. 

To calculate j in from Eq. (|2"6"|) . we have to determine the normalization of the solution 
in physical terms. As a simplification, we consider a complex scalar and a plane wave 
moving around an S 1 parametrized by z G [0, z m ). Going to the rest frame with respect 
to momenta parallel to the brane and reinstating time dependence, we have 

= A/V m(z+t) (27) 

for the plane wave moving towards the UV barrier. To determine the normalization con- 
stant M, we use the standard charge density for a Klein-Gordon particle, j° = Im(0*<9 t 0). 

6 Solutions to Eq. (|23|) in the throat region (z ^> R) are cf) ~ A ^mzJi + 2(mz) + B yJmzYi + 2(mz), 
where A and B follow from the boundary condition on the UV side of the throat and from normalization. 
For sufficiently large z, 4> behaves as <p ~ Acosmz + Bsmmz. The quantization of m is a result of the 
boundary condition for cp or its first derivative at z = z m . 
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It has to be normalized according to 

1 = / dzf J\f = (28) 

JO y mz i~R 

The flux is then given by j- m = j z = lm(<j)*d z <fi). Using the solution of Eq. ( 1271) with the 
normalization of Eq. (j28l . we find 

1 . , 

Jin = — = m m . (29) 

•^IR 



Using this result and Eq. (f24j) . the decay rate of a glueball follows from Eq. ff25|) as 

T ~ m m (m J R) 8+4Z . (30) 

Now that we have the decay rate in the gravity picture, we need to define a vertex 
V in the gauge theory picture which reproduces this result. We model the coupling by a 
term 

Aod 3 Vd^iX-iX^^x^X^Gix) (31) 

in the lOd Lagrangian, where Q denotes the glueball state with canonically normalized 
4d kinetic term. Compactifying the 6 dimensions perpendicular to the brane on a torus 
of size L for the moment and using the KK mode decomposition of Eq. (flOl) . we get the 
effective 4d Lagrangian 

C ^ D E {-\ d ^ **** - \ m * $ * + e 2 ™^ L L Q( X )\ . (32) 

From this, the total decay rate of a glueball into KK modes of the dilaton follows: 

In this formula, pf = pf., is a 4- vector characterizing the momentum of the final-state 
dilaton parallel to the brane, while u>i and ouf are the energies of the initial and final 
state. The 4-momentum of the decaying glueball is denoted by Pi. Introducing the dilaton 
momentum in the compact dimensions as pf ± = 27rn/L, we can replace the sum by an 
integral when we go back to L — ► oo: 

-1 y — (34) 

The decay rate of a glueball into a dilaton is then given by 
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Since the dilaton is massless, ujf = y\pf ± \ 2 + \pf^ | 2 - Going to the rest frame of the glue- 
ball, pi — 0, and performing the momentum integrations, we arrive at 

where we have used Uf = \pf x \ and neglected prefactors of order one. In its rest frame, 
Ui is simply the mass m of the glueball. Comparing with Eq. (1301) . we get 

V ~ y/mjn m 2+2l R 4+21 . (37) 



3.2 Decay rate calculation in the gauge theory picture 

With the effective vertex V at hand, calculating the decay rate of one glueball into gauge 
fields living on a different brane stack is straightforward. Following from Eqs. ( 1321) and 
(1371) . the vertex between a glueball and a KK mode of the dilaton is 

V_ e 2mn { X)/L_ (3g) 

The other vertex in the diagram is still given by Eq. ( TTTT) . Summing over all intermediate 
KK modes, we arrive at an expression very similar to Eq. (fl2j) : 



4+2/ 



y/m m m (mR) v e 



Compared to Eq. (fT2|) . the only difference is the prefactor and the substitution of the 
energy y/s of the colliding gauge bosons by the mass m of the glueball. 

We will analyse Eq. (139]) in two different regimes, namely for m~ l > L and for 
mT 1 L. The former case is the most interesting one from a phenomenological view- 
point. As we argued in Sect. El we can assume that the reheating temperature Trh in 
early cosmology is smaller than L~ x . Accordingly, the mass m of any relic KK modes is 
also restricted by m < L^ 1 . The latter case, on the other hand, can be easily analysed 
in the gravity picture as well. We will perform this cross-check in Sect. 13.31 

For vrC 1 > L, we can make the same simplifications as in Eq. f|T3l) and use Eq. f|T4l) 
for the sum. The decay rate of a glueball into a pair of gauge bosons follows from the 
standard 4d formula: 

r-m^lA^I 2 . (40) 

To get the total decay rate, we have to sum over the iV 2 final state gauge bosons. If we 
denote by Ri and i?2 the AdS scale of the throat containing the initial and the final 
state, respectively, we find 

p8+4« c>8 p8+4Z p8 

r ~ m m m 8+4 < + th_Jk mm m 4 +4 *_ (41) 
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Although the derivation of this decay rate assumed two AdSs x S 5 throats and a torus 
as embedding manifold, it can also be applied to more general geometries, according 
to the discussion in Sect. [2j However, for different throat geometries the dependence 
on the eigenvalues of the angular Laplacian is of course different. These eigenvalues 
entered the discussion through the tunneling probability Eq. (|24p . from which we deter- 
mined the dilaton-glueball vertex in Eq. (1371) . For the example of a Klebanov-Strassler 
throat [1], let us outline how to determine the dilaton-glueball vertex for more general 
throat geometries. Away from the bottom of the throat at r = r s , the warp factor of a 
Klebanov-Strassler throat is well approximated by 

A(r) = 1 + ±!-!L. (42) 

The effective AdS scale R depends on the number of fractional D3-branes at the conifold 
singularity. The metric is still given by Eq. (120]) away from r = r s if one also replaces the 
line element dflf of a sphere by the line element of T 1,1 . For R ^> r ^> r s , which defines 
the throat region, the warp factor is approximately A ~ R 4 ln(r/r s )/r 4 . For r ^> R, 
where the geometry is asymptotically a cone over T ' , we have A ~ 1. Near r = r s , 
the geometry differs considerably from Eqs. (142]) and (1201) and the throat is cut off by 
the Klebanov-Strassler region. For an order of magnitude estimate, one can neglect the 
logarithmic r dependence of the warp factor away from r = r s and approximate the 
Klebanov-Strassler region by a sharp cut off [32]. Thus, the tunneling probability from 
the throat into the conical region can be (approximately) calculated from the effective 
Schrodinger equation, Eq. (123]) . The dependence on the eigenvalues of the Laplacian 
on T 1,1 enters through the potential, where they replace the corresponding eigenvalues 
1(1 + 4) on an S 5 . Moreover, for an AdS warp factor and a sharp cut off, the incoming 
flux is given by Eq. (|29[) . as before. From Eq. (125]) . one can then determine the decay 
rate and match the vertex such that this decay rate is reproduced. 

Let us now consider the case mT x <C L. We will also assume i C L for simplicity. 
Recalling that = 2n\n\/ L and a = 2ttA/L, we can approximate the sum in Eq. (1391) 
by an integral, 



-y 



L 6 f-* m 2 - (2ir) 2 n 2 /L 2 + ie J (2vr) 6 m 2 - p 2 + ie 



(43) 



where p = 2ix n/L. The resulting expression is just the propagator of a massless particle 
in a mixed, energy-configuration-space representation, with the 'energy' m characterizing 
the invariant 4- momentum. This is of course expected in the large L limit, where the 
torus goes over to flat space and the infinite KK tower is replaced by the underlying 
higher- dimensional dilaton field. The integral is evaluated in the Appendix, the outcome 
being 

iA p TTb^ 



(27r) 6 m 2 - p 1 + ie A 2 2K " v ' 

where H^(x) = ^(x) + iY 2 (x) is a Hankel function and we have neglected prefactors 
of order one. Using the asymptotic forms of the Bessel functions for large and small 
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arguments, Eq. (jUJ) can be simplified as follows: 



m2 u+( A^ j^ eimA form- 1 «A 

— Ho{mA)~{ A > (45) 
A 2 yj^ hrm-^A. V ' 

Inserting these results in Eq. fl39l) . we get the matrix elements Ai for these two cases. 
The corresponding partial decay rates follow from Eq. fj40|) . Summing over all final state 
species, we find 

( r x r 2 mmm ii+4i for m -\ ^ A 

1 ~ < R 8+4l R s _ ■ (46) 

[ 1 A8 2 m m m 8+ for m _1 ^> A 

Again, the same discussion as before applies concerning the extension to more realistic 
geometries. As a consistency check, we should examine, whether the appropriate limiting 
cases of Eqs. (l4~T!) and (piS"]) coincide. The regions of validity of the two calculations have 
a common border for A <C mT x ~ L. Indeed, for this choice of parameters the first term 
in Eq. (JUJ) dominates and the result agrees with the second line of Eq. (1461) . 



3.3 Some calculations in the gravity picture 

As in the sections before, we consider two AdSs x S 5 throats embedded in a 6-dimensional 
torus of uniform size L. The geometry is that of a multi-centered black 3-brane, the metric 
being 

ds 2 = A- 1 / 2 (-dt 2 + dx\ + dx\ + dxj) + A 1 / 2 (dx 2 4 + ■■■ + dx 2 9 ) (47) 

with 

= 1 + £ ( ^ + ^ ) ■ (48) 

^{ix^-A.+nL^ \x ± -A 2 + nL\*J 

The positions of the two throats are denoted by A x and A 2 , their AdS scales by R\ and 
i?2- The vector x± refers to the coordinates 24, . . . , xg in the torus. The sum in the warp 
factor A(x±) is due to mirror effects in the torus. Again, this is not a consistent com- 
pactfication. Including O-planes, for example, would give extra contributions to the warp 
factor (see [3]). We try to calculate the transition of a dilaton between different throat 
regions, which is the gravity counterpart to the gauge theory calculation in Sects. 13.11 
and 13.21 The equation of motion for the dilaton is given in Eq. ( !2T|) . Inserting Eq. (T4T1) 
in Eq. (pTj) and using y/g = A 1 / 2 , one gets 

d n d n + A(x ± ) d„ 0" = 0. (49) 

The indices fi and n run from to 3 and from 4 to 9, respectively. Using the 4d Klein- 
Gordon equation, one arrives at 

d n d n (j) + A(x±) m 2 (j) = 0, (50) 

where m is the kinetic energy perpendicular to the branes. Like Eq. ( |23i) . this has the 
form of a Schrodinger equation. Contrary to Eq. (|23j) . however, there is no potential 
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barrier separating the throat region and asymptotically flat space, since the potential 
V = —m 2 A(x±) is strictly negative. The difference comes from using cartesian coor- 
dinates perpendicular to the branes in Eq. (j47j) rather than spherical coordinates in 
Eq. f[2~Ul) . Still, a wave in the throat region, moving away from the horizon, is reflected to 
a large part before entering asymptotically flat space. In cartesian coordinates, however, 
this is due to the steepness of the potential well. 

To determine the transition probability V of a dilaton between two throat regions, 
one has to solve Eq. (150]) with appropriate boundary conditions. Then V is the ratio 
of incoming flux in one throat and outgoing flux in the other throat. In general, the 
corresponding calculation is difficult. However, if the torus is very large (L — > oo) and 
the throats are sufficiently far apart (A ^> m _1 ), the problem effectively splits into two 
simpler calculations. Namely, the latter condition means that the de Broglie wavelength 
of the particle is small compared to the distance of the throats. A transition between two 
throats can then be described as a two-step process. For simplicity, we take the initial 
state in the first throat to be an s-wave. Only a small fraction of the outgoing flux reaches 
the asymptotically flat region, the probability being (cf. Eq. (I24p for / = 0) 

^-(mi^) 8 . (51) 

In between the two throats, one has a free spherical wave, approximating a plane wave 
near the second throat. The absorption cross section (per brane world- volume) for such 
a plane wave was calculated in [21]. Neglecting prefactors of order one, it reads 

ct 2 ~ m 3 Rl (52) 

Near the second throat, the incoming flux will be diluted by a factor of A~ 5 , since the 
free spherical wave is expanding in 6-dimensional flat space. The absorption probability 
by the second throat thus is 

* ~ f ~ «*> 

The transition probability between the two throats is just the product V\Pi- If we denote 
by m m the mass gap in the first throat, using Eqs. ( 1251) and ( 1291) the decay rate from the 
gravity calculation follows as 

r-f^r^W (54) 

This is precisely what we found in Eq. (146 j) for A ^> m~ l and I = 0. The crucial 
ingredient is the A~ 5 dependence. That it agrees in both calculations is, however, not 
too surprising. In the gauge theory calculation, it came from the propagator in a mixed 
energy-configuration-space representation (cf. Eq. (|43p ). The same is of course true in 
the above gravity calculation, although we have not stated it explicitly. 

There is yet another situation where the decay rate between two throats is com- 
paratively easy to obtain. Let us consider only one throat for the moment. We do not 
need to specify the precise form, but will assume that it is finite and reasonably well 
approximated by a slice of AdSs times some compact manifold M.. The prime example 
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certainly is a Klebanov-Strassler throat, whose interpretation as a stabilized Randall- 
Sundrum model was given in [32]. Let us denote by Ri the (approximate) AdS scale 
of the throat and by L the size of the embedding manifold, whose precise geometry is 
again not important. One has L > Ri, since otherwise the throat could not be glued 
into the manifold. If the embedding manifold is of minimal size, L ~ R x , KK modes 
with masses m n <C R^ 1 cannot resolve its precise geometry. We can then describe the 
embedding manifold by the Planck brane in a Randall-Sundrum model. Let us consider 
the Kaluza- Klein expansion of the graviton in the throat. If we restrict ourselves to an 
s-wave with respect to the compact manifold M. multiplying the slice of AdSs, we can 
take the action from [33] obtained in the context of Randall-Sundrum phenomenology 

S = [ d 4 x V | -~d a h$ d a h^ --ml hfrj W v{n) + 4= 7== h$ ) . (55) 
J \ 2 2 VWRi ) 

The effective 5d Planck scale M 5 is determined by Mf ~ Mf Q R\. We have included 
the coupling of the KK modes h!$ to the energy-momentum tensor T^ v on the Planck 
brane, which we will need in a moment. For KK modes with m n <C R[ , the masses are 
determined by 

1' 



Ji(m n /m m ) ~ m n ~ \ n + -j irm m , (56) 

where m m = is the inverse conformal length of the throat (cf. Sect. I3.ip and we have 
used the asymptotic form for large arguments of the Bessel function J\. This is consistent 
for n somewhat larger than 1. The coupling constants g n were calculated in [33], the result 
being 

-((^)) 2 - 1 )" 1/2 ^«" (57) 

In the last step we have used the asymptotic forms for the Bessel function Y\. 

Let us return to the case of two throats and consider another throat in the embedding 
manifold. We take the throat to be AdSsxS 5 such that it can be equally well described 
by a stack of D3-branes. Again, its AdS scale R2 cannot be larger than L, and since we 
have assumed L ~ R±, one has R\ > i?2- The corresponding number N2 of D3-branes 
follows from Eq. (CQ) as N2 ~ Mf Q R\. Now, when viewed from the first throat, the gauge 
theory on the stack of A^ 2 D3-branes resides on the Planck brane. Therefore, the graviton 
KK modes in this throat couple directly to the energy-momentum tensor of the gauge 
theory. Using the last term in Eq. (|55|) . the decay of these KK modes into the other 



7 The usual orbifold boundary conditions were taken for the derivation of coupling strengths and 
masses of graviton KK modes. It is not immediately clear whether the same boundary conditions follow 
from a reduction to 5d of a lOd geometry since the effective theory is defined on an interval instead of 
an S /Z2 orbifold. However, one can rederive the Randall-Sundrum model on an interval if one takes 
Gibbons-Hawking terms [34] at the IR and the UV brane into account. Varying with respect to the 
metric yields a condition similar to the Israel junction condition, to be evaluated only at one side of the 
brane. Inserting the background metric, one finds the relation between the cosmological constants on 
the brane and in the bulk as well as the usual boundary conditions for the fluctuations (see e.g. [35] for 
a derivation of the Israel junction condition using Gibbons-Hawking terms). 
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throat can be calculated as a decay into gauge fields^ By the standard formula, the 
decay rate of a KK mode with mass m n into one species of gauge fields is 



r 




(58) 



There are iVf gauge fields in the adjoint representation of U(A^). Summing and using 
Eqs. ( I57I) and p]), the total decay rate follows: 



This result should be compared with Eq. ( flTi) from the pure gauge theory calculation. 
The distance A between the two throats cannot be smaller than their AdS scales R\ and 
R 2 - Since we have also assumed L ~ Ri and m n <C R^ , the second term in Eq. (|4T|) is 
dominant. Using I = for the s-wave that we have considered and L ~ R\, we get the 
same result as Eq. (1591) . including the factor of m IR ! 

The above process is just the reverse of the energy loss by the heated Planck brane 
considered, e.g., in [27,28]. Our calculation can also be viewed as a rephrasing, using 
partly the gauge theory picture and partly the gravity picture, of the tunneling calcula- 
tion performed in [6] . In these papers, the decay rate of graviton KK modes between two 
throats was calculated in a 5d model with two AdS 5 slices which are glued together at 
a common Planck brane, assuming equal AdS scales R\ = R 2 . However, besides giving 
the corrections due to different AdS scales, from the above derivation it is maybe more 
evident why the result is correct also in a genuine lOd setup. 

The decay rate T ~ (mR) 4 m m from [6] was used in a number of papers [9-12] in 
the context of reheating after brane-antibrane inflation. Moreover, [15] contains a careful 
analysis in a 5d model of effects related to the finite length of realistic throats. In this 
paper, the global KK modes in the two-throat system are determined. Tunneling of KK 
modes is then viewed as the decoherence of wave packets, which are set up in one throat. 
We have used this picture in Sect. 13.11 

Tunneling in a compact lOd setup with throats was considered in [13,14]. For the 
case mT 1 > L, a decay rate of T ~ (mR) w m m was derived, assuming that the particle 
has to tunnel through two barriers described by the potential in Eq. (1231) . We see a con- 
ceptual problem with this approach since we do not know how to justify a 1-dimensional 
quantum-mechanical picture (this 1 dimension being the radial coordinate) in the two- 
throat case. But even if we accept this description for the moment, there are further 
issues related to the two-barriers assumption: The barriers extend to values of r ~ m _1 
as can be seen from Fig. |H Since m" 1 ^> R and r measures the physical distance for 
r ^> R (cf. Eq. ( 12"U|) ). the width of each barrier is given by mT 1 . This just reflects the 
fact that a particle with mass m has a de Broglie wavelength of m _1 . Accordingly, the 
particle has to tunnel through two entire barriers only if the distance A between the two 
throats is ~ 2m~ 1 . Indeed, from Eq. (14TT) for / = and since L > A, we get a decay 

8 There are also decays into the fermions and scalars in the gauge theory. However, the corresponding 
decay rates have the same order of magnitude as the decay rate into gauge fields. 




(59) 
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rate of F ~ (mi?) 16 m IR in this case, in agreement with [13, 14]. However, if A is smaller 
than ~ 2m _1 , the particle has to tunnel through a smaller barrier. Correspondingly, the 
decay rate becomes larger, as can be seen from Eq. (|4T1) . 

The case mT x < L (without assuming m <C L) was also considered in [13,14]. 
It was found that the decay rate can be much larger than T ~ (mi?) 16 m IR if a certain 
resonance condition is fulfilled. We have not determined the decay rate for this case and 
therefore have no result to compare with@ It would be interesting, though, to evaluate 
Eq. (139]) for vrC 1 < L and to see whether one can reproduce the results from [13, 14] as 
well as their resonance condition. 

4 Conclusions and Outlook 

We have determined the energy loss rate p of a throat which is heated to a certain 
temperature as well as the decay rate T of single KK modes localized in a given throat. 
As a simplified setup we have chosen a 6-dimensional torus with two AdSsxS 5 throats. 
However, as we have argued in Sect. [21 our results stay parametrically correct for more 
general embedding manifolds and throat geometries. Especially, they are applicable for 
two Klebanov-Strassler throats if the curvature of the space connecting them is not larger 
than the inverse distance. 

In earlier investigations [6, 9-15] of the decay of KK modes between throats, the 
decay/tunneling rate was determined from solving wave equations in a given gravity 
background. Most results were derived for the simple model of two AdS 5 slices, glued 
together at a common Planck brane. As we have explained in Sect. I3.3[ this calculation 
is difficult to perform in a genuine lOd setup. Inspired by [21], we instead chose the dual 
gauge theory picture for our calculations. Namely, each AdS 5 xS 5 throat can be equally 
well described by a corresponding stack of D3-branes. Both brane stacks are coupled by 
the supergravity fields in the embedding space. The energy transfer rate from a heated 
throat then follows from simple tree-level quantum-field-theory processes. For the decay 
rate of throat-localized KK modes which are dual to glueballs, we first had to determine 
the glueball-supergravity vertex. To this end, we have calculated the decay rate of throat- 
localized KK modes into flat lOd space in the gravity picture. Then, we have determined 
the glueball-supergravity vertex by demanding that the decay rate following from this 
vertex give the same result. We have also presented some cross-checks from the gravity 
picture in Sect. 13. 3[ 

From our analysis, we were able to determine the dependence of the energy transfer 
and decay rates on the distance A between two throats as well as on the size L of the 
embedding manifold. For example, this is relevant for the analysis of reheating after 
brane-antibrane inflation. In such models, one often considers inflation occurring in one 
throat, whereas the standard model branes reside at the bottom of another, longer throat. 
In that way, the generation of the right level of density fluctuations is reconciled with a 

9 Note that we have assumed that L ^> A,m~ x in deriving Eq. (TH)]) . This result is therefore not 
suitable to compare with the results from [13, 14] where the limit of extremely large L was not taken. 
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solution of the gauge hierarchy problem a la Randall- Sundrum. For a viable reheating of 
the standard model sector, it is crucial that the energy from brane-antibrane annihilation 
is transferred efficiently into the standard model throat. This question was analysed 
in [9-12,14]. We find that, as long as the embedding manifold is not of minimal size, 
the energy transfer rate Eq. (|l8p is considerably lower than the rates previously derived 
in [9-12]. Given our results, it will be interesting to reconsider reheating after brane- 
antibrane inflation. 

Our results remain applicable if one deals with a small stack of D3-branes0 An 
interesting setup is the following: Consider that the standard model resides on some 
D-branes in a given Calabi-Yau orientifold. Since they are a common feature of flux 
compactifications, such a manifold will typically contain several throats [5]. Modelling 
the standard model branes by a small stack of D3-branes, we can estimate the rate of 
energy loss to the throats in early cosmology from Eq. (fl8|) . According to Eq. ([1]), with 
iV being small, we just have to replace one factor of R 8 by the corresponding power 
of the lOd Planck scale, M{^ 8 . The throat sectors, which are heated up in that way, 
may provide interesting dark matter candidates [10,14]. Later in cosmological evolution, 
throat-localized KK modes may decay back to the standard model. The correspond- 
ing rate can be estimated from Eq. (T4T1) . again replacing R 8 , by M^ 8 . The decay rate 
strongly depends on the angular momentum of the throat-localized KK modes. We have 
given this dependence explicitly for the angular momentum with respect to an S 5 in an 
AdSsxS 5 throat. Moreover, we have outlined how to determine this dependence for other 
manifolds, e.g. the (approximate) T 1,1 in a Klebanov-Strassler throat. Depending on the 
cosmological epoch, the decaying KK modes may influence the abundances of light el- 
ements or lead to diffuse gamma-ray background radiation, both effects being strongly 
constrained by observations (see e.g. [36]). Along these lines it may even be possible to 
impose certain phenomenological constraints on multi-throat compactifications. 

More generally, one may discuss several cosmological scenarios where reheating takes 
place either in the standard model sector or in a throat (as is the case after brane- 
antibrane inflation) and the standard model resides either at the bottom of a throat or 
somewhere in (the rest of) the Calabi-Yau orientifold. The energy transfer and decay 
rates that we have calculated can then be used in a set of Boltzmann equations to 
determine the evolution of energy densities of the standard model and throat sectors. 
We leave these interesting applications for future work. 

Acknowledgements: We would like to thank J. Braun, F. Brummer, X. Chen, 
D. Dietrich, L. Kofman and M. Trapletti for helpful comments and discussions. 

10 The sole exception is the decay rate of brane-localized states. In this case, our derivation of the 
vertex from the gravity picture does not work, since supergravity is not a good approximation. 



21 



Appendix: Evaluation of the propagator 



In Eq. (|43p . we had to evaluate the following propagator in a mixed, energy-configuration- 
space representation: 



d 6 p 



iA p 



(2tt) 6 m?-p + ie' ^ 

We perform the integral for imaginary values m — *> e m l 2 m and use analytic continuation. 
The integral changes into 



d 6 p 



iA p 



(2vr) 6 m 2 + p 2 ' 

We can then employ the identity c _1 = J °° dre~ CT for Rec > and get 
-1 



(61) 



(2vr) 6 J 
-1 



dr d b p e lAp e 



dr 



dpi e 



[A i pi P -p\'< 



(2?r)6 7 °oo 

-)3 J o T 3 



dp 6 e iMp6 e- p ^ 



(62) 



(47T) 



We have used that A 2 = A\ + • • • + A\. According to Eq. 3.471.9 in [37], this integral 
can be evaluated in terms of the modified Bessel function K_ 2 = K 2 , which yields 



—1 m 2 
(2vr) 3 A 2 " 



K 2 {mA). 



(63) 



Following from Eq. 9.6.4 in [38], K 2 is related to the Hankel function H 2 = J 2 + iY 2 . 
The above expression can be written as 



% m 



(4tt) 2 A 2 



H+(e in/2 mA). 



(64) 



The Hankel function has a branch cut along the negative real axis. Therefore, one can 
analytically continue back to real values m — *> e~ ln ^ 2 m, which gives 



-i m 



(4vr) 2 A 2 



H+(mA). 



(65) 
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